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The random magnetic flux problem on a lattice and in a quasi one-dimensional (wire) geometry is 
studied both analytically and numerically. The first two moments of the conductance are obtained 
analytically. Numerical simulations for the average and variance of the conductance agree with the 
theory. We find that the center of the band e = plays a special role. Away from e = 0, transport 
properties are those of a disordered quantum wire in the standard unitary symmetry class. At the 
band center e — 0, the dependence on the wire length of the conductance departs from the standard 
unitary symmetry class and is governed by a new universality class, the chiral unitary symmetry 
class. The most remarkable property of this new universality class is the existence of an even-odd 
effect in the localized regime: Exponential decay of the average conductance for an even number of 
channels is replaced by algebraic decay for an odd number of channels. 
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I. INTRODUCTION 

The. concepts of scaling0~i and of the renormalization 
groupQ have provided crucial insights into the localiza- 
tion properties of a. quantum particle in a random but 
static environments Beyond a typical length scale de- 
pending on the microscopic details of the disorder, the 
localization problem can be described by an effective field 
theory that is uniquely specified by the dimensionality of 
space and theJundamental symmetries of the microscopic 
Hamiltonian.B Correspondingly, the disorder is said to 
belong to the orthogonal, unitary and symplectic ensem- 
bles, depending on whether time reversal symmetry and 
spin-orbit coupling are present or not.uu 

However, not all disordered systems belong to one of 
these three standard symmetry classes. One example is 
the Integer Quantum Hall Effect, for which the scaling 
theory in the unitary universality class cannot explain the 
observed jumps in the Hall resistance,!!] since it predicts 
that all states are localized in two-dimensions. Instead, 
a new scaling theory was proposed for the Integer Quan- 
tum Hall Effect, where, in addition to the longitudinal 
conductivity that controls the scaling flow in the uni- 
tary ensem ble , .the Hall conductivity appears as a second 
parameter .liBO 

In this paper we consider a different example. It is 
the so-called random flux model, which describes the 
localization properties of a particle moving in a plane 
perpendicular to a static magnetic field of random am- 
plitude and vanishing mean.Ej lj In the literature, dif- 
ferent points of view have been offered with regard to 
the localization properties and the appropriate symme- 
try class of the random flux problem. In Refs. p"8|-|24| 
it has been claimed that, since the magnetic field has 
a vanishing mean, the only effect of the random mag- 



netic field is to break time reversal invariance, and hence 
that the localization properties are those of the standard 
unitarv-symmetry class. On the other hand, Zhang and 
ArovasEZI have argued that this argument might be too 
naive and that a scaling theory closely related to that 
of the Kosterlitz-Thouless transition controls the local- 
ization properties of the random magnetic flux problem. 
They predicted that states are localized in the tails of 
the spectrum whereas close to the center of the band a 
line of critical points of the Kosterlitz-Thouless type is 
formed. Related point of views can be found in Refs. 
p5|-|32|. Finally, it has been proposed in Ref. ^ that the 
random flux model shows critical behavior at the band 
center e = only, whereas its localization properties are 
those of the unitary ensemble for energies e ^ 0. 

In the third scenario, the behavior at e = is gov- 
erned by an additional symmetry, the so-called chiral or 
particle-hole symmetry. The chiral symmetry can also 
be found in the related problem of a particle hopping 
on a lattice with random (real) hopping amplitudes.^ 
In the one-dimensional version of this problem, it is well 
established that the ensemble-averaged density of states 
diverges at the band center e = OEjEZI and that the en- 
semble averaged conductance decays algebraically with 
the length L of the systemO For comparison, in the uni- 
tary symmetry class, the density of states is continuous 
at e = 0j23 while the conductance decays exponentially 
with L. (The one-dimensional random- hopping problem 
has been studied in many incarnations, cf. Refs. [hi] |46[ ) 
For two-dimensional systems, the effect of,the chiral sym- 
metry was studied by Gade and WegnercZl (see also Refs. 
[jq-pq) . They argued that the presence of the chiral sym- 
metry results in three new symmetry classes, called chiral 
orthogonal, chiral unitary, and chiral symplectic. For dis- 
ordered systems with chiral unitary symmetry, all states 
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are localized except at the singular energy e = at which 
the average density of states diverges. The relevance of 
the chiral unitary symmetry class to the random flux 
problem was pointed out by Miller and WangH (Only 
the chiral unitary class is of relevance, since time-reversal 
symmetry is broken in the random flux model.) 

For the two-dimensional random-flux problem, suffi- 
ciently accurate numerical data are notoriously hard to 
obtain. Although a consensus has emerged that states 
are localized in the tails of the spectrum, it is impossi- 
ble to decide solely on the basis of numerical simulations 
whether states are truly delocalized upon approaching 
the center of the band, or only deceptively appear so as 
the localization length is much larger than the system 
sizes that are accessible to the current computers. More- 
over, it is very easy to smear out a diverging density 
of states in a numerical simulation (compare Refs. [l^j22] 
and [34] ). In short, no conclusion has been reached in 
the debate about the localization properties of the two- 
dimensional random flux problem. 

Here, we focus on the simpler problem of the random 
flux problem on a lattice and in a quasi one-dimensional 
geometry of a (thick) quantum wire with weak disorder, 
and restrict our attention to transport properties, no- 
tably the conductance g. For a wire geometry, numerical 
simulations can be performed with very high accuracy, 
and very good statistics can be obtained. Moreover, pre- 
cise theoretical predictions for the transport properties 
can be made, both for the unitary symmetry class, and 
for the chiral unitary symmetry class. The wire geometry 
allows us to quantitatively compare the analytical predic- 
tions for the various symmetry classes and the numerical 
simulations for the random flux model. This compari- 
son shows that, away from the critical energy e = 0, the 
L-dependence of the average and variance of the conduc- 
tance g are those of the unitary ensemble. At the band 
center e = 0, (g) and varg are given by the chiral uni- 
tary ensemble. Hence, we unambiguously show that in a 
quasi one-dimensional geometry, the localization proper- 
ties of the random flux model are described by the third 
scenario above, in which the e = is a special point, 
governed by a separate symmetry class. Although our 
theory is limited to a quasi one-dimensional geometry, it 
does show the importance of the chiral symmetry at the 
band center e = and may thus contribute to the de- 
bate about the localization properties of the random flux 
problem in higher spatial dimensions. 

This paper was motivated by twOr*ecent works. First, 
in a recent paper, one of the authorsEa computed (g) and 
var g numerically for the random flux model in a wire 
geometry to a very high accuracy. While for nonzero en- 
ergies e, the result was found to agree with-analytical cal- 
culations for the unitary symmetry classj5H~Ell for e = a 
clear difference with the unitary symmetry class was ob- 
served. Second, for the chiral symmetry classes, a scaling 
equation for the distribution of the transmission eigenval- 
ues in a quasi one-dimensional geometry was derived and 
solved exactly in the chiral unitary case by Simons, Al- 



tland, and two of the authors .E3 This scaling equation 
is the chiral analogue of the so-called-J)orokhov-Mello- 
Pereyra-Kumar (DMPK) equation,E5L3 which describes 
the three standard symmetry classes and was solved ex-. 
actly in the unitary case by Beenakker and Rejaei.Ea 
However, for the chiral unitary case, analytical results 
for the L dependence of (g) and var g were lacking, so 
that a comparison between the theory and the numerical 
results of Ref. [34| was not possible. In the present work 
this gap is bridged. 

In a wire geometry, the chiral unitary universality 
class undergoes a stpkiag even-odd effect first noticed 
by Miller and Wang£3c3 The conductance g decays ex- 
ponentially with the length L if the number of channels 
N is even, while critical behavior is shown if N is odd, 
even in the limit of large N that we consider here. In 
the latter case, the average conductance (g) decays alge- 
braically, while the conductance fluctuations are larger 
than the mean. We analyze how the even-odd effect fol- 
lows from the exact solution of the Fokker-Planck equa- 
tion of Ref. [32] and compare with numerical simulations 
of the random flux model. 

We close the introduction by pointing out that the ran- 
dom flux problem is also relevant to some strongly corre- 
lated electronic systems. In both the Quantum HalL-Efk 
feet at half-fillingLaLa and high T c superconductivityj!3'L3 
strong electronic correlations can be implemented by 
auxiliary gauge fields. In this context, the random flux 
problem captures the contributions from the static trans- 
verse gauge fields. Notice that the chiral symmetry is not 
required on physical grounds both for the Quantum Hall 
Effect at half-filling and for high T c superconductivity. 
Another area of applicability, for .our results is the pas- 
sive advection of a sc^a&AsilaOEj'Lj and non-Hermitean 
quantum mechanics. OT3Eal£3 Finally, the striking sen- 
sitivity of the localization properties in the random flux 
problem to the parity of the number N of channels is 
remarkably similar to that of the low energy sector of a 
single antiferromagnetic spin- N/2 chain to the parity of 
A,l3 on the one hand, or to the sensitivity of the low 
energy sector of N coupled antiferromagnetic spin-1/2 
chains to the parity of NfLU on the other hand. 

The paper is organized as follows. The random flux 
problem in a wire geometry is defined in section [n| The 
average and variance of the conductance are calculated 
analytically in section III. Analytical predictions are 



compared to the numerical simulations in section IV . We 
conclude in Sec. 



II. THE RANDOM MAGNETIC FLUX MODEL 

In the random flux model one considers a spinless elec- 
tron on a rectangular lattice in the presence of a random 
magnetic field with vanishing mean. The magnetic field 
is perpendicular to the plane in which the electron moves. 
In this paper, we study the random flux model in a wire 



2 



(a) 








e , 2 


e ll2 


e 2 , 2 


e 3 , 2 


e 4 , 2 









e ,i 


9i,i 


e 2 ,i 


e 3 ,i 


e 4a 





.7 = 3 

3 = 2 
3 = 1 



m = 



m = 



(b) 







y C R 








c i < 




, r R 


v,- 







L 



FIG. 1. (a) Lattice with TV = 3 threaded by random 
magnetic fluxes Q m ,j in the disordered region < m < M. 
(b) Quantum wire with a disordered region of length L = Ma. 



Incoming amplitudes are c„ + and 
amplitudes are _ and c„/ 



1,. 



whereas outgoing 
, N , in the left and 



right leads, respectively. In a quasi one-dimensional geometry, 
M > N. 



geometry and for weak disorder. This system is described 
by the Hamiltonian 



Hip 



-t[lp m +l,j + 1pm-l,j] 

-t(l- S ]tN )e iem - 1 ip mt3+1 
-ta-^iJe-^-^mj-i, 



(2.1) 



where ipm,j is the wavefunction at the lattice site (m,j), 
labeled by the chain index j = 1, . . . , N and by the col- 
umn index m, see Fig. 0(a). The Peierls phases 9 m j re- 
sult from the flux Q m ,j = 9m+i,j — 9 m ,j through the pla- 
quette between the sites (m,j), (m+l,j), (m + l,j + l), 
and (m, j + 1). (The flux O mj - does not uniquely de- 
termine all the phases along all the bonds. We have 
used this freedom to choose the nonzero phases along 
the transverse bonds only.) 



We consider a system with Hamiltonian (2.1) where 
the phases & m ,j take random values in a disordered strip 
< m < M only, and are zero outside.E£l We assume 
that the disordered region is quasi one-dimensional, i.e., 
M N 3> 1, corresponding to a thick quantum wire. In 
the disordered region, the Peierls phases 9 m j are cho- 
sen at random in such a way that the magnetic flux 
mj = 0m+i,j—9 m ,j is uniformly distributed in [— pir,pir] 
with < p < 1. To be precise, with 9 m j given, 9 m +i.j 
is chosen from the interval [0 m ,j — P^, 9 m ,j + P^} with 
uniform probability l/2pw. The parameter p controls 
the strength of disorder. We assume weak disorder, i.e., 
p«l. 



The boundary conditions in the tra nsve rse directions 
that are implied by the Hamiltonian (2.1) are "open", 
i.e., there are no bonds between the chains j = 1 and 
j = N. In this case, Ti has a special discrete symmetry, 
called the particle-hole or chiral symmetry: Under the 



transformation ip„ 



(-lf'+^mj, one has H -> -H. 



Hence, for each realization of the random magnetic flux, 
the chiral symmetry ensures that there exists an eigen- 
state of H with energy — e for each eigenstate of H with 
energy +e. Note that the band center e — is a special 
point. The chiral symmetry is broken by the a ddit ion of 
a random on-site potential to the Hamiltonian (2.1). An- 
other way to break the chiral symmetry is to add bonds 
between the chains j — 1 and j = N and to impose peri- 
odic boundary conditions in the transverse direction for 
iV odd. The presence of the chiral symmetry may have 
dramatic consequences for charge transport through the 
disordered wire, as we shall see in more detail in the next 
sections. 

In order to find the conductance g of the disordered re- 
gion with the random flux, we first compute the transfer 
matrix A4. To the left and to the right of the disordered 
region, the wavefunction ^> m ,j that solves the Schrodinger 
equation Ti%p = sip can be written as a sum of plane waves 
moving to the right (+) and to the left (— ), 



i/=l ± 
v=l ± 



sin— r, m < 0, 



sin k u "" N + 1 

,±jfc„m ^ l ,j n 

sin , m > M. 

smk„ N + l 



where cos/c^ = —e/2t — cos[istt/(N + 1)]. The prefactor 
1 / sin k v is chosen such that an equal current is carried 
in each channel. The number N c is the total number 
of propagating channels at the energy e, i.e., the total 
number of real wavevectors k v . We are interested in the 
transport properties for e close to 0, where N c = N, 
and ignore the distinction between N c and N henceforth. 
The coefficients c^+and c^ ± arc related by the transfer 
matrix M [see Fig. ||(b)], 



(2.2) 



Note that M. v y is a 2 x 2 matrix in Eq. ([L4 ). Current 
conservation requires 



(2.3) 



where S3 = 03 ® Iat, CT3 being the Pauli matrix and 
In the N x N unit matrix. In addition, at the special 
point e = 0, the chiral symmetry of the Hamiltonian (2.1) 
results in the additional symmetry 



SiTWSr =M, 
where Ei = o\ ® Ijy. 



(2.4) 
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The eigenvalues of MM 1 * , which occur in inverse pairs 
exp(±2xj), determine the transmission eigenvalues Tj = 
1 / cosh 2 Xj and hence the dimcnsionless conductance g 
through the Landauer formulaHEa 



JV 



5 = E^ = E 



cosh 



(2.5) 



In the absence of disorder, all exponents Xj are zero, and 
conduction is perfect, g = N. On the other hand, trans- 
mission is exponentially suppressed if all Xj s tire larger 
than unity. The smallest Xj determines the localization 
properties of the quantum wire. 

For the quasi one-dimensional geometry M ^> N ^> 1 
that we consider here and on length scales much larger 
than the mean free path associated to the random mag- 
netic field, the microscopic details of the microscopic 
Hamiltonian TL should no longer be important. Rather, 
the crucial ingredients are the symmetries of TL. For 
nonzero energy, the only sym metry of M. is given by 
current conservation, Eq. (2.3). In this case, for quasi 
one-dimensional systems with sufficiently weak disorder, 
the probability distribution P(xi, . . . ,xjv;L) of the pa- 
rameters Xj is governed by the so-ca!kid_porokhov-Mello- 
Pereyra-Kumar (DMPK) equation,Elt^ 



dL 4N ^ fl"-- 



dxj 



(2.6a) 



J = Y[ | sinh 2 Xj - sinh 2 x k \ 2 | sinh(2a;j)|. (2.6b) 

k>j k 

Here L = Ma is the length of the disordered region, a 
being the lattice constant. The mean free path I depends 
on the disorder strength and on the deta ils of the micro- 
scopic model. The derivation of Eq. ( [2.6D assumes i 3> A, 
A being the wave length at the Fermi energy. The initial 
condition corresponding to perfect transmission at L = 
is P{x\ , . . . , xn; 0) = Yl 3 S(xj ) . The Fokker-Planck equa- 
tion ( |2.6D describes the unitary symmetry class. For 
e = 0, in add ition to current conservation, the chiral 
symmetry (2.4) has to be taken into account. In Ref. |62| 
it was shown that for weak disorder (p <C 1) the distri- 
bution P(xi, . . . , xn', L) s atisfies again a Fokker-Planck 
equation of the form (2.6), but with a different Jacobian 



,dP 
dL 



1 N 

— Y 

2N ^ 



3=1 



d_ 

dxj 



j = n \ 8inii ( x j - x k)\ 2 - 

k>j 



(2.7a) 
(2.7b) 



This Jacobian describes the chiral unitary symmetry 
class. As was shown in Ref. and as we shall see in 
more detail in the next se ction , as a result of t he rc - 
placement of the Jacobian ( 2.6b| ) by the Jacobian ( 2.7b| ), 



the statistical distribution and the L-dependence of the 
conductance g at energy e — is quantitatively and qual- 
itatively different from that away from e — 0. In Ref. |62| 
it was shown that there exists a quantum critical point 
induced by the randomness when TV is odd within the 
chiral unitary symmetry class. Away from zero energy, 
the transport properties of the disordered wire are those 
expected from the standard unitary symmetry class. A 
derivation of Eq. (2.7) is given in Appe ndix 

The physical picture underlying Eqs. ( |2.6| ) and (2.7) is 
that the parameters Xj undergo a "Brownian motion" as 
the length L of the disordered region is increased. The 
Jacobian J describes the "interaction" between the pa- 
rameters Xj in this Brownian motion process. The key 
difference between the unitary case and the chiral uni- 
tary case is the presence of an inter action with "mirror 
ima ged" e igenvalues Xj in Eq. ( 2.6b ), which is absent in 
Eq. (2.7b). To see this, we note that both for the unitary 
and for the chiral unitary cases, the Jacobian J vanishes 
if a parameter Xj coi ncid es with x k , k ^ j. However, 



in the unitary case (2.6b), J also vanishes if Xj coin- 



cides with a mirror image — x k , k ^ j, or if Xj = (i.e., 
Xj coincides with its own mirror image). The vanishing 
of the Jacobian J implies a repulsion of the parameters 
Xj in the underlying Brownian motion process. Hence, 
whereas Xj feels a repulsion from the other N — 1 param- 

; feels 



eters x k , k j , in the chiral unitary case ( |2.7| ) 
an additional repulsion from the N — 1 mirror images 
—Xk, k 7^ j, and from its own mirror image — Xj in the 



standard unitary case, 
It can be shownEjl 



that the parameters Xj repel 
each other by a constant force in the large-L limit, ir- 
respective of their separation. This long-range repulsion 
results in the so-called "crystalization of transmission 
eigenvalues" : The fluctuations of the parameters x j are 
much smaller than the spacings between their average 
positionsEJ Away from zero energy, i.e., in the unitary 
symmetry class, all Xj can be chosen positive because 
of repulsion fspm their mirror images, and their average 
positions are£3 



( Xj ) = (2j-1)L/2N£, j = l,...,N. 



(2. 



In the chiral unitary symmetry class, the Xj can be both 
positive and negative since there is no repulsion from the 
mirror images, and one hasE2l 



{ Xj ) = (N + l-2j)L/Nt, j = l, 



(2.9) 



In the unitary symmetry class and in the chiral unitary 
class with even N the net force on each parameter xj is fi- 
nite, and they grow linearly with the length L. Hence, by 
Eq. ( |2.5| ), the conductance g is exponentially suppressed 
for L N£. However, for the chiral disordered wire 
with an odd number of channels N, the net force on 
the middle eigenvalue a?(jv+i)/2 vanishes: it remains in 
the vicinity of the origip_and the conductance is not ex- 
ponentially suppressed.E3 Thus, the quantum wire with 
random flux with an odd number N of channels goes 
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through a quantum critical point at zero energy whereas 
it remains non-critical for an even number N of channels. 
A more quantitative description of this even-odd effect is 
developed in the next section. 



III. MOMENTS OF THE CONDUCTANCE 

A. Method of bi-orthonormal functions 

To calculate the moments of the conductance g, we 
make use oLthe exact solution of the Fokkcr-Planck equa- 
tion Q,B 



P(xx,...,x N ;L) oc JJe 



WZ *i 



x II ( x j ~ x k) sinh(o;j - x k ). (3.1) 



The proportionality constant is fixed by normaliza tion of 
the probability distribution. A derivation of Eq. (3.1) is 
presented in Appendix [b]. 

The moments of O-Can be computed from the n-point 
correlation functionsEl 



Rn{x\ : • ■ ■ j X n: -L) 

Nl f + °° 

W^nyJ-oo dXn+1 



(3.2) 



dx N P(xi, . . .,x N ;L), 



and the Landauer formula (2J5). For example, the first 
and second moments of g are 



(fj) 



dx 



Ri(x;L) 



x cosh 2 x 



, ... , , , , R 2 (x 1 ,x 2 ;L) 
{g ) = j dxi j dx 2 l2 



' cosh xi cosh X2 



+°° R.jx-L) 
dx - r -. 

-oc cosh X 



(3.3a) 



(3.3b) 



Here we compute R n (xi, . . . , x n ; L) using the method of 
bi-orthonormal functions developed by MuttalibEa and 
FrahmEII for a disordered wire in the unitary symmetry 
class. The idea is to construct, for any given N and L, a 
function Ki,(x,y) with the following properties, 



dx Ki,(x, x) = N, 



(3.4a) 



-Hoc 



dyK L {x,y)K L {y,z)=K L {x,z), (3.4b) 

P({xi}\L) = c N det [KL(xi,Xj)] iJ=1 ^ N . (3.4c) 

If such-a function exists, it is known from random matrix 
theoryEl that c N = l/N\ and 



R n ({xi};L) = det [K L (x l ,x 



Our construction of the function Kl{x, y) st arts with a 
representation of P(x\, . . . , xn', L) in Eq. (3.1) as a prod- 
uct of two determinants. Making use of the identities 



J\(xk - Xj) = det 

j<k 

j J sinh(a;fc — Xj) = det 



J- 1 



j,k=l,...,N 



lJN+l-2j)x k 
2 e 



j,k=l,...,N 



we find 



P({xi};L) cx det [</>j( x k)] jik=1 , 

x det {xk)] j^k=i,...,N 



where 



4>j (%) 
r)j(x) 



(3.6a) 

(3.6b) 
(3.6c) 



Note that the way we write P as a product of two deter- 
minants in Eq. ( |3.6[ ) is not unique. In particular, we are 
free to replace the sets of functions {<f>j} and {rjj} by an 
arbitrary set of linear combinations {<j>j} and {%}. This 
freedom is crucial for the construction of the function 
Kl(x,v), as we shall see below. 

Since the product of two determinants equals the de- 
terminant of the product of the corresponding matrices 
and since transposition of a matrix leaves the determi- 
nant unchanged, it is tempting to identify K^(x,y) with 



Y^,f=i ^j{ x )Vj(v)- I n tlu s way, Eq. ( |3.4c] ) is satisfied. 
Howe ver, w ith th is choice, the remaining two conditions 
(|3.4aj ) and (3.4b) are not obeyed. This problem can be 
solved by making use of the above-mentioned freedom 
to replace the sets of functions {4>j} and {rjj} by linear 
combinations {4>j} and {fjj}. One easily verifies that if 
we choose these-4inear combinations such that they are 
bi-orthonormalH3 



dx4>j(x)fi k (x) = Sjk, j,k=l,...,N, (3.7) 



all three conditions (3.4) arc met if we set 



3 = 1 



(3.8) 



The construction of the bi-orthonormal functions (f>j and 
rjj is done below. 

First, we define the set {fjj(x)}, j = 1, . . . , iV, by com- 
pleting the square in the exponent of rjj (x) and then nor- 
malizing rjj (x) , 



fjj(x) 



2-kcj 



~{x-£ja) 2 /2a 



(3.9) 



3 ' H,j— 1,. . . ,n 



(3.5) where we abbreviated 
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a = L/N£, 



£j=N + l- 2j. 



(3.10) 



B. Average conductance 



The functions <fij , being linear combinations of <pj (x) = 
x 3 ~ , are polynomials themselves, too. Their (maximal) 
degree is N — 1. As a first step towards their construction, 
we define the polynomials 



dy(iy/a) j - 1 e-( y+ix ) /2<T , (3.11) 
which satisfy the special property 

dx Pj (x)fj k (x) = {e k y-\ (3.12) 



Notic e th at Pj(x) is of degree j — 1. According to 
Eq. ( 3.12 ), the overlap matrix between the polynomi- 
als pj and the Gaussians f)j is independent of L. Con- 
struction of bi-orthonormal functions <fij and fjj is thus 
achieved by choosing L-independent linear combinations 
of th e po lynomials pj that diagonalize the overlap ma- 
trix ( 3.12 ). Shis is done using the Lagrange interpolation 
polynomials^ 



L m {x) = Yl 



(3.13) 



which are of degree N — 1 and obey L m (e n ) — 6 mj7l . We 
infer that the desired polynomials (pj (x) are given by 



4>j(x) 



1 

2na 



dyL {iy/a)e~ [y+lx)2/2 ' 7 . (3.14) 



Putting everything together, we find that 

N J. + OQ 

Kl(x,z) = J dyLjiiy/v) 



3=1 

x exp 



(y + ix) 2 + (z — Sjcr) 2 
2a 



(3.15) 



Now, moments of the conductance g can be calculated 
with the help of Eq. ( |3.5| ). In particular, we find that the 
average and variance of g are given by 



(9) 



var g 



-oc 



dx 



K L (x,x) 



(3.16) 



cosh x 

hoo r+°° K L (x 2 ,x 1 )K L (xi,x 2 ) 
axi / ax2 

30 J — OO 

+ °^^. (3.17, 
-co cosh X 



cosh 2 x\ cosh 2 X2 



After some shifts of integration variables and with the 
help of the Fourier transform of cosh -2 x, 



dx 



cosh 2 x sinh(7ry/2) 



k=l 



(3.18) 



we obtain from Eqs. ( |3.15| ) and ( 3.16| ) an expression for 
the average conductance (g) at the energy e = that 
involves one integration and one (finite) summation only, 



N 



(9) = E 



7/2 



(3.19) 



rn — 1 



dtj 



L m (e m -iy)ye- ( -y+^) 2 °/2 
2sinh(vry/2) 



where, as before, a — L/Ni. In the lim it N ^> 1 at fixed 
a (the so-called thick-wire limit), Eq. (|3.19D can be fur- 
ther simplified. Hereto we use the second identity of Eq. 
( |3.18| ) to cancel the Lagrange interpolation polynomial 
in the coefficient c m , 



dy _ 



y 2 a/2 



A, 



3 71" 

{-m + 1, 



n 



feeA,, 



ly - 



2k 



-m + N}. 



(3.20) 



In the limit N — > oo, only m's close to (JV + l)/2 con- 
tribute to (<?). For those m, w e may replace the infinite 
product on the r.h.s. of Eq. (|3.20D by unity, and find 
c.,„ = (2/ira) 1/2 . Hence, for JV > 1 even, 

(g) = jA^o^a/Tr) = Jl. V e -™ V / 2 , (3.21a) 

V 7TCT V 7T(T * — ' 



whereas for N 3> 1 odd, 



(ff> = \/ — ^ 3 (0|2ia/7r) 

7T<7 



-7T/T ^ ' 



-m 2 a/2 



(3.21b) 



Here $2 and $3 are the Jacobi's theta functions^ 

The dramatic difference between even and odd chan- 
nel numbers discovered in Refs. |3^,^ follows immediately 
from Eqs. (|3.19|) or ( |3.2l| ) in the regime L 3> N£. For 
even N, each term in the summation decays exponen- 
tially with L. The exponential decay of (g) is governed 
by th e slowest decaying terms in the summation in Eq. 
(|3.19[ ), i.e., the contributions from e m — ±1, i.e., from 
to = N/2 or to = N/2 + 1. Hence for L > N£ we find 



(9) 




TV even. (3.22) 
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FIG. 2. Average conductance at zero energy as a function 
of L/£ (£ = AT£) for a quantum wire with a finite number 
of channels N = 1, 2, 3, 4 and in the chiral unitary symmetry 
class. The even-odd effect is clearly visible for L/£ ^> 1: (i) 
Exponential decay of (g) for N even, (ii) Algebraic decay of 
(g) for N odd. 
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FIG. 3. Average conductance as a function of L/£ for a 
quantum wire in the limit N 3> 1. Curves (a) and (b) are for 
large odd and even N at zero energy, when the system is in 
the chiral unitary symmetry class. Curve (c) shows (g) for a 
quantum wire in the unitary symmetry class with the same 
mean free path as in (a) and (b). 



Eq. ( 3,^2 ) allows us to identify £ as the localization 
lcngt h£3 F or odd N, there is one term in the summa- 
tion ( 3.19 ) that does not decay exponentially with L. 
It is the contribution from the channel with e rn = 0, 
m = (N + l)/2. In this case, we again define £ = Ni, 
although it is now merely a crossover length scale, to 
be the characteristic length scale above which the slow 
algebraic decay of (g) sets in, 




N odd, 



(3.23) 



In Fig. g we have shown the average conductance for 
]V = 1,2,3,4 as a function of L/£ and the asymptotic 
result for large N. 

To study the average conductance in the diffusive 
regime £ <C L <C £, we use the Poisson summation for- 
mula 



^5(o;-2m-l) = ^e i ^- 
Tiez 

£5(x-m) = 5> 2 ™ 



1) 



nez 

jinx 



in6Z 



to convert Eq. (3.21) into 

(.9) 



L 



(3.24a) 
(3.24b) 

(3.25) 



¥ £ (-l)"e _7r " i/2L , N > 1 even, 



N > 1 odd. 



Hence the even-odd effect is non-perturbative in L/£ and 
we see that £ = N£ is the characteristic length scale at 



which the even-odd effect shows up. Whereas the leading 
terms are identical, the first non-perturbative correction 
to (g) differs by a sign for even and odd N. 

In Fig. H we plot the average conductance for N ^> 1 as 
a function of L/£ and compare with the unitary symme- 
try class, which is appropriate for energies e ^ 0. InJthe. 
unitary symmetry class, (g) decays exponentially ,E3~Eil 
irrespective of the parity of N, but with a different local- 
ization length £ u , 



(g) cx e 



£u = 2N£. 



(3.26) 



The unitary symmetry class is appropriate for the ran- 
dom flux model if the energy e is nonzero. Hence, moving 
the energy e away from zero causes a factor 2 increase in 
the localization length if the number of channels is even, 
and a dramatic decrease in the average conductance if N 
is odd. 



C. Variance of the conductance 

Proceeding as in the previous subsection, we find from 
Eqs. pl7|), and fOl), 



JV 



varg = > <■„,.,, ( 

m,n— 1 
N 



<j/2 



(3.27) 



with the coefficients 



V — C 



Vl L m (e n -i yi )e-^+^ 2 ^ 2 
2sinh(7ryi/2) 
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FIG. 4. Same as Fig. g, but for the variance of the con- 
ductance. 



FIG. 5. Same as Fig. |3|, but for the variance of the con- 
ductance. 



dy 2 



dy 



2sinh(7n/ 2 /2) 
(y 3 + Ay) L m (e m - ly) e -(v+*.»>V2 



12sinh(7ry/2) 

We plot varg, which is computed from Eq. ( [3. 27 ) for 
N = 1,2,3,4, in Fig. ||, together with the thick wire 
limit N 3> 1. The even-odd effect is clearly seen when 

In the limit N — > oo at fixed N£/L further simplifica- 
tions are possible. We find 



var g 



^ , ^ t fm.nfn.m + ^ f f m 
T 



m— — oo n— — oo 



m— — oo 



f 




-rrTcr/2 



(to — n) erf (rn-\/ <r/2j — \m — n\ 



18lT<T 



™2 | -1\ -mV/2 

to + a e ' , 



(3.28) 



where the primed summations are restricted to even 
(odd) to and n for N odd (even). The error function 
erf(x) is defined as 



2 f x 

erf (x) = —= I die . 
V 71 " Jo 



For L > £ Eq. fl3.28|) simplifies to 

^e- L /«, TV »1 even, 

AT > 1 odd. 



var g : 



(3.29) 



The variance of the conductance decays exponentially for 
large even N with the same decay length as the aver- 
age (g), while varg decays algebraically for large odd N. 
Note that (g) and varg decay with the same power of L. 



After some tedious algebra starting from Eq. (3.2S) 
to extract an expression well suited for an asymptotic 
expansion in small L, we find for the diffusive regime 



2_ 

15 



var g 



15 3 \L J 



N 3> 1 even , 

A" > I odd. 

(3.30) 



Again, we see that the difference between even and 
odd channel numbers shows up in terms that are non- 
perturbative in L/£. The leading term 2/15 in varg is 
universal and twice the value of its counterpart for _a. 
disordered quantum wire in the unitary ensemble. LjEil 
Hence, moving the energy e away from zero decreases 
the conductance fluctuations by a factor two in the dif- 
fusive regime. The factor two decrease of var g upon 
breaking the chiral symmetry is reminiscent of the fac- 
tor two difference for the conductance fluctuations be- 
tween the standard orthogonal and unitary symmetry 
classesS The enhancement of the conductance fluctua- 
tions at e = had been observed previously in numerical 
simulations of the-two-dimensional random flux problem 
by Ohtsuki et al.EZI 

Figure |^ contains a plot of varg versus arui-oi. 
fers a comparison with the unitary symmetry class.EHtil 
In the unitary symmetry class, varg takes the univer- 
sal value 1/15 in the diffusive regime L <C N£, while 
varg cx exp(— L/2£ u ), £ u = 2N£, in the localized regime 
L->N£. 



IV. NUMERICAL SIMULATIONS 

In this section we present numerical simulations for the 
conductance g in the random flux model (2.1). The aver- 
age and variance of g were studied previously by Avishai 
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et alM and by Ohtsuki et alM for the random flux model 
in a square geometry. However, for a comparison with the 
theory of Sec. [II and to identify the symmetry class it 



is necessary to study a wire geometry and large system 
sizes. This is done below. 

For each disorder configuration, we calc ulate the con- 
ductance using the Landauer formula (2.5), which we use 
in the more conventional form 



E 



(4.1) 



Here _/V c is the number of propagating channels in the 
leads and t is the N c x N c transmission matrix, which re- 
lates the amplitudes of the incoming and outgoing waves 
on the left and the right of the disordered sample. The 
eigen values Tj of the matrix tt^ are the same as in Eq. 
(2.5). (The simulations are aimed at energies e close to 
zero, where N c equals N. Hence, as before, we drop the 
notational distinction between N c and N). 

The transmission matrix t is cpmeuted through the 
recursive Green function method.EilE3 In this method, 
N x N matrix Green functions Fj k for reflection and 
Gjk for transmission through the disordered region are 
computed using the recursive rule, 



F(m + 1) = [s - H m - i 2 ^(to)]" 1 , 
G(m + 1) = -tG(m) F(m + 1), 



(4.2) 



where the matrix elements of H m are 



H 



3,k 



j+i,k 



-ta-^ije-*-'- 1 ^--!,*- 



The initial conditions at m = are those of a Green 
function at the edge of an isolated perfect lead: 



F jk (0) = G jk (Q) 
2 



N ■ 



N 
v=l 



e ik " sin 



vjn . vkir 
■ sm ■ 



N + l N + l' 



(4.3) 



where cos k u = —e/2t — cos[^7r/ (N + 1)], see Sec. [n[ The 
scattering channels are those modes with real wavevec- 
tors k v . The Green function that we need is obtained by 
taking into account the perfect lead boundary condition 
on the right of the disordered region, 

F(M) = [FiO)- 1 - t 2 F(M - 1)] _1 , 
G(M) = -t G(M - 1) F[M). 

The matrix Green function G(M) describes the propaga- 
tion from ui — to xtx = M. The absolute value of the 
transmission matrix element at energy e is then given by 



tn,u\ — 4 sin fc M sin k v 
2 



x 



N 



N ■ 



1 ^ 

j,k=l 



Gjk(M) sin 



N + l N + l 



This procedure is repeated for each disorder configura- 
tion, and the average and variance of the conductance are 
obtained by taking an average over 2 x 10 4 samples. The 
transverse boundary conditions are those of Eq. (|2.l|), 
i.e., open boundaries, unless explicitly indicated other- 
wise. We present the numerical results as a function of 
L/£, where £ i s the characteristic length entering Eq. 
flS.22 ) and Eq. ( 3.25 ). We determine £ by compari ng the 
num erical data for L ^> £ to the asymptotic results (3.22) 
and ( ]3~23|V 

Figure |6| shows the average and variance o f th e con- 
ductance at e = of the random flux model (2.1) with 
iV = 15 and N — 16 and disorder strength p = 0.2. When 
L ^S> £, (g) decreases algebraically for N = 15 whereas 
it decays exponentially for N — 16. This is precisely the 
even-odd effectSt 2 ] that we discussed at length in the last 
section. We find excellent agree men t between the numer- 
ical data and the theory of Sec. Ill, which is indicated by 



the solid (odd N) and dashed (even N) lines in the figure. 
The characteristic length £ that gover ns the crossover to 
the slow algebraic decay of Eq. ( 3.23 ) is estimated to be 
280a for N = 15. The localizati on le ngth £ that governs 
the exponential decay of Eq. ( |3.22 ) is estimated to be 
283a for N = 16. As in the case of the average conduc- 
tance, for varg, the even-odd effect can be clearly seen 
for L > £, where the numerical data coi ncide with the an- 
alytic result in the large N limit, Eq. (3.2?). The slight 
discrepancy at very small L happens at M ~ N and may 
be understood as a crossover from quasi one-dimensional 
to quasi two-dimensional behavior. This type of one- 
to two-dimensional crossover was reported in a numeri- 
cal work by Tamura and Ando.cZl The Fokker-Planck ap- 
proach employed in sections || and III is specifically de- 
vised for a quasi one-dimensional geometry and is there- 
fore inapplicable to describe the regime M < N. 

In Figs. and || we consider the dependence of (g) and 
var g on N, p, and e. Figure 0(a) shows the average and 
variance of g for odd N at e = for two choices of N and 
p. We see that the numerical data show fairly good agree- 
ment with the analytic large N result (solid lines) for the 
three cases we examined. For larger disorder strength p, 
the deviations from the analytical result (3.28) is more 
prominent, the stronger disorder data being closer to the 
onset of quasi two-dimensional behavior for small L. The 
agreement between the numerical data for p = 1 and the 
theory of Sec. Ill for L ^ £ is remarkable, in view of 



the fact that the theory was derived under the assump- 
tion of weak disorder, whereas p = 1 corresponds to the 
strongest possible disorder in the random flux model. 

Figur e |7| (b) shows (g) and varg for the random flux 
model (ETj) at e = for an even number N of chan- 
nels. We show the data for three cases (N,p) = (16, 1), 
(16,0.2), and (32,0.2). In the last example we used peri- 
odic boundary conditions in the transverse direction in- 
stead of the open boundary conditions of Eq. fl2.ip . Since 
N is even, the periodic boundary conditions do not de- 
stroy the chiral symmetry, so that the system remains in 
the chiral unitary symmetry class. We see that the re- 
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FIG. 6. The average (a) and the varia nce (b) of the con- 
ductance g for the random flux model (2.1) at e = for 
TV = 15 (circle) and TV = 16 (diamond) with disorder strength 
p = 0.2. For these parameters, we find that £ = 280a for 
TV = 15 and £ = 283a for TV = 16. The solid (dashed) lines in 
(a) are the theoretical result (3. It) for (g) for TV = 15 (16); 
the solid (dashed) line s in (b) are the large odd (even) TV 
analytical results (3.25) for vaxg. 
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FIG. 7. The average and the variance of the conductance 
g for the random flux model (2.1) at e = for odd TV (a) 
and even TV (b). (a) The circles, squares, and triangles are 
the results for (TV,p) = (15,1.0), (15,0.2), and (31,1.0), re- 
spectively. The characteristic length £ is numerically found 
to be £ = 23.7a, 280a, and 39.8a, respe ctive ly. The solid lines 
are the large odd TV analytical result ( |3.28] ). (b) The circles, 
squares, and triangles are the results for (N,p) — (16,1.0), 
(16,0.2), and (32,0.2), respectively. In the last case, periodic 
boundary conditions (p.b.c.) in the transverse directions are 
used, whereas for the first two cases open boundary condi- 
tions are assumed. The localization lengths £ are found to be 
£ = 27.2a, 283a, and 476a, respectively. The white solid lines 
are the large even TV analytical result for the chiral-unitary 
class of section til. 
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FIG. 8. Average and variance of g for disorder strength 
p — 1 away from the critical energy e = 0. The circles and 
squares are the data for TV = 15 and e = 0.02T (f u = 47.8a) 
and for N = 16 and e = 0.02t (£ u = 50.6a), respectively. The 
diamonds and triangles are the data for TV = 15 and e — 
(f„ = 45.0a) and for TV = 32 and e = O.lt (£ u = 1041a), which 
are calculated for the periodic boundary condition. For odd 
TV, periodic boundary conditions break the chiral symmetry. 
The white solid lines are the large TV analytical result for the 
unitary class, taken from Refs. |59|-|6l[ 



suits of numerical simulations are indistinguishable from 
the theoretical curves (solid lines) for both (g) and var g 
except in the quasi two-dimensional regime M < TV. We 
conclude from Fig. ^|(a) and ^|(b) that the localization 
properties of the random flux model at e — are gov- 
erned by the chiral unitary universality class, indepen- 
dent of the disorder strength. 

In figure || we show some results where the chiral sym- 
metry is broken. In this case, charge transport is no 
longer governed by the Fokker-Planck equation (2.7) for 
the chiral unitary symmetry class, but by the Fokker- 
Planck equation (|]^) that is valid for the standard uni- 
tary class.E^TEJ In the figure, numerical results are shown 
for three cases away from the critical energy as well as 
for one case where the chiral symmetry does not exist 
because of the periodic boundary condition imposed for 
odd TV. With the exception of very short lengths, where 
the system becomes quasi two-dimensional, all the data 
for (g) and for var g agre&,with the theoretical prediction 
for the unitary class EBtEI The results indicate that the 
small nonzero energy e — 0.02t is sufficient to cause a 
change from the chiral unitary symmetry class at e = 
to the standard unitary symmetry class. Another inter- 
esting feature to note is that the localization length £ in 
Fig. Sis roughly twice as large as that in the chiral case 
(Fig. 0). (For example, compare the two cases e = 
and e = 0.02t for TV = 16 and p = 1, where we find 
£ = 27.2a and £ u = 50.6a, respectively.) This behav- 
ior was observed earlier in Refs. |88|,[33f[54|. This result 
is consistent with the analytic result that £ differs by a 



factor of 2 between the chiral (£ = N£) and the unitary 
class (£ u = 2N£), assuming that the mean free path de- 
termined by the short-distance physics is identical in the 
two classes. (For the numerical results we may expect 
that the mean free path should not have strong energy 
dependence on the scale of |e| < O.lt.) 



V. CONCLUSION 

In this paper we studied transport properties of a par- 
ticle on a rectangular lattice in the presence of uncorre- 
cted random fluxes of vanishing mean. This problem is 
commonly known as the random flux problem. We con- 
sidered a wire geometry and weak disorder and showed 
that the symmetries of the random flux problem have 
dramatic consequences on the statistical distribution of 
the conductance g. If the energy e is away from the band 
center e — 0, the system belongs to the standard unitary 
symmetry class, while at e = 0, transport is governed 
by an additional symmetry of the random flux model, 
the particle-hole or chiral symmetry. We have compared 
numerical simulations of the average and variance of the 
conductance g in the random flux model in a thick quan- 
tum wire to analytical calculations for the standard uni- 
tary and the chiral unitary symmetry classes, and found 
good agreement for e ^ and e = 0, respectively. 

There are important differences between the conduc- 
tance distribution in the chiral unitary symmetry class 
and the standard unitary symmetry class, both in the 
diffusive and the localized regime. These differences are 
summarized in Table ||. The most striking feature of-the. 
chiral unitary symmetry class is the even-odd effect&tS 
If the number of channels TV in the wire is even, the av- 
erage conductance (g) decays exponentially with length 
L in the localized regime L ^> Nt, whereas for odd TV, 
the decay of (g) is algebraic. The sensitivity to the chi- 
ral symmetry in transport properties is very strong. For 
example, removing the chiral symmetry by a change in 
boundary condition is sufficient to change the universal- 
ity class to the standard unitary one, even in the thick- 
wire limit TV 3> 1. 

Although our theory is limited to a quasi one- 
dimensional geometry and cannot account for the 
crossover from one to two dimensions, it does show 
the importance of the chiral symmetry to the trans- 
port properties of the random flux problem. Taking 
the prominent role played by symmetry for the random 
flux model in quasi one-dimension as a guideline, we 
speculate that a similar picture is appropriate for the 
two-dimensional random flux problem. Hence, following 
Gade and Wegncr^Zl and Miller and Wanga we expect 
that the localization properties of the two-dimensional 
random flux problem are controlled by the unitary sym- 
metry class away from the band center e = 0, whereas 
the band center e = plays the role of a critical energy. 
The random flux problem would thus share with the Intc- 
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TABLE I. Average and variance of the conductance g of a 
disordered quantum wire with N ^> 1 channels and mean free 
path £ in the unitary and chiral unitary symmetry classes, 
for the diffusive regime L <C Nl and for the localized regime 
L 3> Nl. The results for the unitary ensemble are taken from 
Refs. E<l-|6l|. 



ger Quantum Hall Effect, and with the problem of Dirac 
fermions in a random vector potential the existence of 
a single critical energy that lies between energies with 
localized states. There are however two important differ- 
ences with the Integer Quantum Hall Effect. First, there 
is no symmetry that fixes the value of the critical energy 
in the Integer Quantum Hall Effect, while the chiral sym- 
metry of the random flux model implies that criticality 
occurs at the band center e = 0. Second, in contrast to 
the smooth density of states in the Integer Quantum Hall 
Effect, one expects that the density of states in the ran- 
dom flux problem is singular at e = 0. Such a singularity 
of the density of states at e — was observed in the sin- 
gle chain random hopping problem, is suggested by the 
numerical simulations of Rcfs. [49|]34| , and is consistent 
with Gade's analysis of the twttdimensional non-lincar- 
cr model with chiral symmetryjjj and with exact results 
on the pir©|biem of Dirac fermions in a random vector 
potential. E3£il (The latter problem shares the same chi- 
ral symmetry as the random flux problem although it 
differs from the random flux problem in that the mag- 
netic fluxes are strongly correlated on all length scales.) 
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APPENDIX A: DERIVATION OF THE 
FOKKER-PLANCK EQUATION 

In this paper, we described the transport properties of 
a quantum wire in the chiral unitary symmetry class in 
terms of its transfer matrix M. . Our theoretical analysis 
was f ocu sed on a solution of the Fokker-Planck equa- 
tion (2/7) that governs the //-evolution of the probabil- 
ity distribution P{x\, ... ,xn;L) of the eigenvalues e ±2xj 
of AiAi^ . A derivation of this Fokker-Planck equation 
from a different microscopic model was presented in Rcf . 
|62] , Here we present an alternative derivation of Eq. (2.7) 
that is closer in spirit to derivations of the Fokker-Planck 
equation for the unitary class existing in the liter ature.cil 

For the statistical distribution of the parameters Xj, 
the symmetries of the transfer matrix M.{e) are of fun- 
damental importance. For the random flux model, there 
are two symmetries (cf. Sec. ||): 



M(s) S 3 M\e) = E 3 , flux conservation, (Al) 
EiA^(e)Ei = M(-e), chiral symmetry. (A2) 



Here the transfer matrix M. is defined in Eq. ( |2.2| ) and 
Ej = <T,j ® Ijv, where Oj is the Pauli matrix (J = 1,3) 
and 1 at is the N x N unit matrix. 

Because,©! flux conservation (Al) Ai(e) can be param- 
eterized asc3 



M 



Mn Ml2 

M 2 i M22 

U \ /coshA sinhA 
W ) ^sinhAT cosh A 



V 
V 



(A3) 



where U, W , V, and V" are N x N unitary matrices and 
A is a diagonal matrix containing the parameters Xj on 
the diagonal. We are interested in the case of zero en- 



ergy, when the chiral symmetry (A2) results in the fur- 
ther constraints U —W and V = V. Notice that in this 



case, with the parameterization (A3), the parameters 



are uniquely determined by Ai. This is an important 
difference with the unitary symmetry class, where each 
xj is only defined up to a sign. As a result, in the uni- 
tary class, the distribution P(xi, . . . , xjv; L) has to be 
symmetric under a transformation xj — > —Xj for each j 
individually, while no such-symmetry requirement exists 
in the chiral unitary classO 

As the length L of the disordered region is increased 
(see Fig. ^|), the parameters Xj, j — 1, . . . , N are sub- 
jected to a Brownian motion process: As L is increased 
by an amount SL, the parameters Xj will undergo a (ran- 
dom) shift Xj — > Xj + Sxj . We first seek the appropri- 
ate Langevin equations that describe the statistical dis- 
tribution of the increments Sxj. Hereto we note that 
the transfer matrix M. = A4(0;L + SL) is the product 
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SL 




L 



FIG. 9. A thin slice of length SL with a < SL < t < L 
is added to the disordered region of length L. 

of the individual transfer matrices M. = A4\0;L) and 
M 1 = A4(Q; SL) for wires of length L and 5L, respec- 
tively: 



M =MM'. 
We also use that the matrix 

2M n M\ 2 = U sinh(2X)W t 



(A4) 



(A5) 



is hermitian and has eigenvalues sinh2xj, j = l,...,N. 
Hence we find that 

2M u M{ 2 = U(smh2X + 2A)U\ 

A = W f (M n M\ 2 - M u M\ 2 ) U. (A6) 

Making use of the symmetries of M! and of the parame- 
terization ( [A3] ) we can rewrite A as 

A = coshXVM'^M'^V 1 sinh AT 
+ sinhA:V7W' 12 7W' 1 t 2 V t coshX 
+ coshA:V7W' 11 Al' 1 t 2 V t cosh AT 
+ sinh X VM' 12 M'^ sinhX. (A7) 

We take the length SL of the added slice small compared 
to the mean free path £. Within the thin slice the disorder 
is assumed to be uncorrelated beyond a length scale of 
the order of the lattice spacing a <C 5L. In this case one 
has M' = 1 + 0{5L) 1 / 2 , so that the matrix A is of order 
{SL) 1 / 2 itself and we can treat it in perturbation theory. 
As a result, we find that the addition of the slice of width 
SL results in the change 

sinh 2xi - sinh 2xj = 2A„ + 4 V A jkA kj 

J J JJ ^ sinh 2x, ;- sinh 2x k 



or equivalently 



Sxi 



A, 



+ 0(<5L 3/2 ), 



A 2 j sinh 2xj 



(A8) 



cosh 2xi 



cosh 2x4 



+ 2 



AjfcAfej 



0(SL 3 / 2 ) 



(sinh 2xj — sinh 2x k ) cosh 2xj 



(A9) 



It remains to find the first two moments of Aj k . Hereto 
we make an ansatz for the distribution of the transfer 
matrix M.' . Because A4' is close to 1, it is natural to 
parameterize it in terms of its generator, 



M' = exp A. 



(A10) 



From the symmetry requirements ( |A1| ) and (A2) we de- 
duce that A has the form 



A=iV <g>l 2 + W ®ai, 



(All) 



where V and W are hermitian NxN matrices. We choose 
a convenient statistical distribution of M! by assuming 
that V and W have independent, Gaussian distributions 
with zero mean and with variance 



(VijVki) = (WijWki) = Mjk^. 



(A12) 



Then we find that the first two moments of A are given 
by 

SL 

(Ajfc) = 5jkShah(2xj)—, 
SL 

(A jk A k:j ) = cosh 2 ^ + x k)jjj- 



Combining this with Eq. ( |A9| ), we conclude that under 
addition of a narrow slice of width SL -C £, the parame- 
ters Xj undergo a shift Xj — > Xj + Sxj with 



{Sxj) SL = j-r p J2 C0th ^ 
k^j 
SL_ 



Xk), 



(SXjSx k )sL = -rr^jk, 



(Al3a) 



(A13b) 



all higher moments vanishing to first order in SL. Equa- 
tion (A13) is equivalent to the Fokker-Planck equation 



APPENDIX B: SOLUTION TO THE 
FOKKER-PLANCK EQUATION 

In this appendix, we p rese nt an exact solution for the 
Fokker-Planck equation (2.7), closely following the exact 
solution of the DMPK equationaa the unitary symmetry 
class by Beenakker and Rejaei.E3 We start by rewriting 
Eq. as 



dL 2N ^ dxi 

3=1 J 



dP 

dxj 



2P 



dx-j 



(Bla) 
(Bib) 



3<k 



where the initial condition is 
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A' 



(Blc) 



P(xi, ...,x N ;0) = Y[ S(xj). 



The key step towards the exact solution of Eq. jBl] ) is 
the transformation 



p({x j y,L) = 



j^J sinh(xj 

j<k 



Xfc) 



*({*,-}; L), (B2) 



which changes the Fokker-Planck equation (Bl) into a 
Schrodinger equation, 



dL ~ 



2N^ dx] + O.N 



1 N 
— Y 



a 2 * 

2N ^ ~dxj 



2N 



an 

dxi 



d 2 n 

dxf 



(B3) 



Here U = (N - l)(N - 2)/6 + (TV — l)/2. Thus, 
\&(xi, ■ ■ ■ , a: at; L) obeys a Schrodinger equation in imagi- 
nary time L that describes N identical free particles on 
the line, — oo < x < oo. (For comparison, in the unitary 
symmetry class, one finds that \& obeys a Schrodinger 
equation for N identical particles moving in the presence 
of a potential oc sinh - 2x which repels the x's away from 
the origin £3) 

Since the probability distribution P(x\, . . . ,xn',L) is 
symmetric un der a permutation of the Xj's, it follows 
from Eq. (B2) that \&(xi, . . . , Xn; L) must be antisym- 
metric, i.e., it must describe the imaginary-time evolution 
of N identical fermions. At L — 0, the initial condition 
(Blc) implies that all Xj coin cide at the origin. Hence, 
at L = 0, the transformation (|B2) is singular. We-avoid 



this problem by starting with the initial condition 

1 N 

p({xjh o|{2/ fc i) = m II 5 ^ - y*u))> 

CT j = l 



(B4) 



where all the initial values are different, and send e to 
zero at the end of the calculation. The summation is 
over all permut ation s er of 1, . . . , N. 

To solve Eq. (B3), we denote by G(x;L\y) the single- 
particle Green function of the diffusion equation obeying 



? dG _ 1 d 2 G 
~dL~2N Ih? 1 



G(x;0\y)=6(x-y). (B5) 



Solution of Eq. ( p5\ ) yields 



G(x;L\y) = 



2ttL 

Then the Slater determinant 
1 



e -m(-~v) 2 



(B6) 



*({*,-}; L|{j/ fc }) = M det [G(x f ,L\y k )}.^ h 



x e 



-ULj 



(B7) 



is antisym met ric in xi, . . . , xn and obeys the Schrodinger 
equa tion (B2). Using the inverse of the transformation 
(B2), we obtain that 



P({x j };L\{y k }) = 
n{*jY,L\{Vk})U 



sinh(xj — Xfc) 
sinh(yj - y k ) 



(B8) 



is the solution to the Fokker-Planck equation (Bl) with 
the regularized initial condition (B4). 

We finally take the limit e — > 0. This limit must be 
treated with care in view of the denominator of Eq. ( |Bg| ) . 
With the help of 



det 



jvf r :2 

2L 



Nie 
~2L 



j,k=l,...,N 

N(N-l) 



(B9) 



j<k 



Xk 



0(e 2 



the singularity oc e N ( N coming from the denomi- 
nator in Eq. (B8) is cancelled. We thus recover Eq. (|3.1|). 
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